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Abstract 

We discuss the determination of polarized parton distributions from 
a NLO analysis of recent experimental data. We extract the first 
moment of the polarized quark and gluon distribution and assess the 
corresponding uncertainties, with special regard to those related to 
scheme dependence. 
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We discuss the determination of polarized parton distributions from a NLO 
analysis of recent experimental data. We extract the first moment of the 
polarized quark and gluon distribution and assess the corresponding un- 
certainties, with special regard to those related to scheme dependence. 



The theoretical and experimental information on inclusive polarized deep- 
inelastic scattering has considerably improved recently. On the one hand, the 
polarized structure function gi has now been measured over a reasonably wide 
range in x at several scales, both for proton and deuterium targets. ^^"^ On 
the other hand, a complete determination of the two-loop polarized anomalous 
dimensions is now available. ""^'^ This makes a fully-fledged NLO analysis of 
polarized parton distributions both feasible and justified. Such an analysis is 
necessary in order to extract the physically relevant matrix elements of quark 
and gluon operators from the measured gi{x, Q"^). 

The simplest leading-twist polarized observables are the nucleon matrix 
elements of the quark axial currents ^ = for ^ target of mass M 

and spin the matrix elements {s\j^j\s) = Mai{Q'^)s^ are related to the first 
moment ofTi{Q'^) = J^dx gi{x,Q'^) by 

ri(O') = i Eriie-aa. = 5 [(e')Cs(0')ao(Q') + CNs(Q')aNs(0')] , (1) 

where the average quark electric charge is (e^) = '^^7=1^1^ ^^'^ singlet 
and nonsinglct axial charges arc respectively oq = Y^^i=i ~ 
Ym=\ i^i ~ (^^)) ^i' ^^"^ C'(Q^) H-rc perturbatively computable coefficient fimc- 
tions, which at NLO are given by Cs(g2) = Cns(Q^) = l~^ + 0{al).'^ While 
the nonsinglet axial charge aNs(Q^) only depends on scale due to heavy quark 
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thresholds and is otherwise determined by SU(3) symmetry in terms of hy- 
peron /3 decay constants, the singlet charge measures a nontrivial combination 
of the first moments of the polarized quark and gluon distributions: 

Cs(Q2)ao(Q2) = C,(l, g')AE(l, Q^) + Cg{l, Q^)Ag{l, Q^) (2) 

where AS(A^, Q^) = Yll=i j^dx x^~^[l^qi + A^i](x, Q^) is the moment-space 
polarized singlet quark distribution, other Mellin moments are defined simi- 
larly, C,(1,Q2) = i + 0{a,{Q'')) and (7^(1, Q^) = 0(a,(Q2)). 

A direct determination of ao(Q^) would be possible only if gi were known 
for all X at one scale; since this is not the case we can only determine it 
indirectly through eq. (2) from quark and gluon distributions, the primary 
quantities which enter the QCD evolution equations. All the available ex- 
perimental information on gi{x,Q^) is then summarized by constructing a 
set of polarized parton distributions (gluon and quark singlet and nonsin- 
glet) at a reference scale Qq, from which gi can be obtained for all x and 
by solving the evolution equations. We parametrize these distributions as 
Af{x, Qo) = Nfr]fx"f {l—x)^f (1+o/a;) (with Nf chosen in such a way that the 
first moment of A/ equals rjf). We evolve these distributions to all the points 
in the {x, Q^) plane where 51 (x, Q^) has been measured, and thereby determine 
the best fit values of the parameters (further details of the fitting procedure 
are given in refs. 7, 8). We use the NLO form of anomalous dimensions and 
coefficient functions, determined in the MS rcnormalization scheme. 

The choice of factorization scheme is particularly subtle in the polarized 
case because of the extra ambiguity related to the prescription used to define 
the 75 matrix in dimensional regularization; without such a prescription the MS 
scheme is not fully defined. We will use the so-called Adler-Bardeen scheme,^ 
which is defined in terms of a specific factorization scheme change on the singlet 
anomalous dimensions and coefficient functions of refs. 4,5 (which arc in turn 
obtained from MS factorization and the 't Hooft-Veltman prcsc;ription for 75 
by a scheme change which enforces the conservation of the nonsinglet axial cur- 
rent). The AB scheme is constructed by imposing conditions on first moments 
of singlet anomalous dimensions and coefficient functions, whose factorization 
scheme dependence is especially delicate. 

To understand this, recall that oq depends on because the singlet axial 
current, unlike its nonsinglet counterpart, is not conserved in the chiral limit 
in the quantized theory due to the axial anomaly: S^js = «/fjtre^'^'"^F^yFp<j. 
As a consequence, the current is multiplicatively renormalized: = 
^^^\oia)j^, with an anomalous dimension 7^^^ = ^^^^(x^g + ^(^s) that starts 
at two loops. The anomaly equation then implies that the (gluonic) anomaly 
operator mixes with the axial current according to ■^Uf^tre'^'^^'^Ffi^Fpa- = 
J^^Has)dfj_j^ . At the level of operator matrix elements, the current is 



2 



parametrized by ao(Q^), and the partonic quantity which mixes with it can be 
identified with the first moment of the polarized gluon distribution: 



(3) 



This is possible because at LO asAg{l, Q"^) is scale independent. 

Eq. (3) imphes that the scale dependence of Ag(l,Q^) causes a 
strong scheme dependence of AE(1,(5^): upon a general scheme change 



(Tg') = [1 ^nd since a,Ag{l,Q^) is scale in- 



dependent at LO it follows that AS' — AS does not vanish asymptotically 

if Zqg 0. Because at LO AS is scale-independent this implies that even 
asymptotically the ambiguity on AS(1, Q^) is of the same order as its size, i.e. 
the definition of AS(1,(5^) is entirely ambiguous. This is not the case for the 
gluon: even if we redefine it by setting Zgg ^ then Ag'(l,(5^) — Ag(l,Q^) 
is asymptotically constant, but since in the same limit Ag(\^Q^) diverges (as 
aj^) the relative scheme ambiguity in Agf(l,(5^) vanishes asymptotically. 

The definition of AS(1,Q^) is thus intrinsically arbitrary and can only 
be fixed on the basis of physical requirements. A simple choice is based on 
the observation that since the quantities on either side of eq. (3) evolve in 
the same way, their difference is scale independent: we can then identify this 
scale-independent combination with the singlet quark first moment 



to all orders. This choice has several advantages: a) only when AS(1) is scale 
independent does the Ellis- JafFe ansatz, which would equate it to the scale- 
independent octet axial charge as, aquire an unambiguous meaning beyond 
LO;^ b) the polarized quark distribution can then be defined directly in terms 
of a hard observable (such as a jet cross section), by properly factorizing all 
soft contributions to the cross section into the parton distributions;^*^ c) AS(1) 
eq. (4) coincides with the matrix element of the canonically defined quark 
helicity operator. 

The AB scheme choice exploits the fact that eq. (3) is automatically satis- 
fied with the factorization choice of refs. 4,5, and enforces the definition eq. (4) 
of the polarized quark distribution. This still leaves the freedom to further re- 
define the normalization of A(/(l,(5^) (i.e. to vary Zgg in the scheme change), 
which can be fixed entirely at NLO by requiring that eq. (3) be also true at 
NNLO: in fact, the Adler-Bardeen theorem^^ implies that there exists a scheme 
where eq. (3) is valid to all orders. Knowledge^'^ of the 0(0;^) contribution to 
7^^) then fixes entirely the scheme choice for the first moments of parton distri- 
butions. As a consequence, in the AB scheme the evolution of ao, as given by 
eq. (3), is automatically correct through NNLO when the parton distributions 



AS(l)=ao(02)+n/f^A5(l,Q2) 



(4) 
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are evolved at NLO. The factorization scheme for all other moments is finally 
fixed by insisting that the full scheme change from the scheme of refs. 4,5 to 
the AB scheme is independent of x. 

The growth of the first moment of the polarized gluon distribution with 
scale implies sizable scaling violations for moments close to the first, which 
evolve at LO by mixing with the gluon distribution. Measuring these scaling 
violations then determines the polarized gluon distribution. This determina- 
tion is only asymptotically scheme independent, while at finite it is subject 
to a factorization scheme ambiguity.* On the other hand, large scaling viola- 
tions driven by a large gluon distribution also imply large scheme ambiguities, 
in particular in the determination of a^. For instance, it is clear from eq. (2) 
that at NLO ao(Q^) will be affected by a scheme ambiguity of order OfgAg, 
which only decreases with scale as Ug. 

Fitting to the data of refs. 1,2 we find® 

ao(10 GeV^) = 0.14 ± 0.10 (exp.) (th.) 
AE(1) = 0.5±0.1, A5(l,10GeV2) = 3.0± 1.6. 

The largest source of experimental error comes from the uncertainty in the 
behaviour of the partons in the unmeasured small x region, ^'^'^^ which could of 
course in principle lead to an infinite (either positive or negative) contribution 
to first moments. In rcf.8 it is constrained by our specific choice of fimctional 
form for the initial parton distributions: assuming that the small x exponent 
a f can be fitted from the data amounts to assuming that the asymptotic small 
x behaviour has already started setting in in the measured region {x > 0.003). 
The data turn out to determine reasonably well the small x behaviour of the 
nonsinglet {a^s = —0.68 ±0.15), but only very poorly that of the singlet quark 
(as = 0.41 ± 0.38) and gluon (ctg = -0.47 ± 0.30); this is then reflected in a 
sizable statistical uncertainty in the first moments. The bulk of the theoretical 
error is due to lack of knowledge of higher order corrections, as manifested 
by scheme dependence, which we estimate by varying the renormalization and 
factorization scales by a factor of two around Q^. Further sources of error. 



Note that the AltareUi-Parisi equations allow a determination of all moments of Ag from 
scaling violations, including the first, provided only that the input distributions are sufficiently 
smooth. This is despite the fact that the singlet first moment evolves multiphcatively, and thus 
measurements at different scales determine only one linear combination of AS(1) and Ag{l). 
However, two functions of x whose moments are all identical, but whose first moments are 
different, can only differ by a term proportional to S{x). Such terms are not experimentally 
measurable, even in principle, and are thus by construction not included in the definition 
of parton distributions. The first moment Ag(f) is thus uniquely determined by analytic 
continuation from Ag{N), / 1. In practice this analytic continuation is performed implicitly 
through the choice of the form of the parametrization of the initial parton distributions. 
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which include the variation of the position of heavy quark thresholds, the 
value of as and SU(2) and SU(3) breaking (as reflected in the values of the 
nosinglet axial charges) are negligible in comparison. 

Since the error on the determination of the polarized gluon distribution is 
dominated by the statistics, it can be improved when more precise or abundant 
data in the presently measured region become available: for instance, includ- 
ing the data of ref. 3 in the analysis leads to A(7(l, 10 GeV^) = 2.6 ± 1.3, even 
though the statistical weight of these data is relatively modest. The determina- 
tion of the axial charge oo, on the other hand, is dominated by the systematics 
related to higher order corrections, and the uncertainty in the small x extrap- 
olation. The former could be improved by higher statistics data which provide 
a more complete coverage of the x — plane (allowing evaluations of the first 
moment at essentially a single value of Q^)- However, to obtain a significant 
improvement in the latter, data in the unexplored small x region will be re- 
quired. A substantial improvement in the determination of the singlet axial 
charge will thus only be accomplished by the next generation of experiments, 
such as would be possible at HERA with a polarized proton beam.^^ 
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